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A direct design method for solving the problem of 
robustness to cross-coupling perturbations in multivariable 
control systems is presented. The method uses numerical 
optimization procedures to manipulate the system feedback 
gains as direct design variables. The manipulation is accom- 
plished in a manner that produces desired performance by 
pole placement and robustness by modification of the minimum 
singular values of the system return difference matrix. 

Channels affected by cross-coupling perturbation may be 
recognized by the character of their transfer function plot. 
The mechanism used by the pole placement and robustness 
routine in obtaining a ro bust - design is evident from the 
gain changes associated with the transfer function diagram 
and the zero shifts shown on pole-zero plots. The pole 
placement and robustness routine uses gain equalization and 
zero assignment to modify the . characteristics of the system 
in the areas of low singular values, producing a robust 
design. 

A modification of the pole placement and robustness 
routine that may be applied to the design of robust 
observers is also presented. Using feedback and filter gains 
as direct design variables a practical design procedure for 
robustness recovery in observer based systems is obtained. 
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I. INTRODUCTION 

With the rising interest in multivariable control theory 
brought on by increasingly complex systems the need has 
arisen to develop design methods that will allow the 
designer tc specify system performance while at the same 
time ensuring relatively high stability margins or robust- 
ness. In the single-input single-output (SISO| case the 
designer has had the tools to do these tradeoffs in the form 
of Nyguist, Bode and root locus plots. In the multi-input 
multi-output (MIMO) case the classical methods are not 
totally appropriate. 

With the increased interest in MIMO systems numerous 
methods of design have been employed to obtain suitable 
system performance and robustness with varying degrees of 
success. One primary method of design is to keep the plant 
as decoupled as possible throughout the design so that each 
individual element may be controlled independently and 
designed essentially as a single loop system. Rosenbrock 
[Ref. 1 ] has developed a procedure where the multilccp 
system is modified into a system that has diagonal 
elements that are much larger than any off-diagonal 
elements. This diagonally dominant system is then in a form 
where conventional Nyguist type techniques can be employed 
in the analysis. A third common MIMO design method is that 
of the linear Quadratic (L Q) method. This method uses a 
quadratic cost functional and optimization principles to 
allow the designer tc design for various performance levels 
by adjusting the matrix weighting terms used in the cost 
function. The major difficulty with all of the above 
methods is that they are not necessarily robust . This is 
especially true for cross-coupling terms between loops. 
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The primary achievement of this thesis has been the 
incorporation of the time domain pole placement design 
procedure with a method of using return difference matrix 
singular values to improve the robustness of the design. The 
technique, which utilizes a modern optimization routine, can 
significantly assist the designer in obtaining robustness in 
the face of cross-coupling perturbations. It has also been 
shown that the cross-coupling perturbation problem can be 
detected by using classical open-loop Bode diagrams as well 
as modern control analysis. The pole placement and robust- 
ness design routine developed for this thesis has been used 
cn several problems discussed in recent literature. In these 
studies the pole placement and robustness design code has 
proven capable of meeting the desired goals cf pole 
placement and robustness and also brought to light seme 
interesting aspects of the cross-coupling perturbation 
problem. A slightly modified pole placement and robustness 
routine has proven effective in the design of robust 
observers. 

The remainder of the thesis will present background 
material on SISO systems in Chapter Two and on MIMO systems 
in Chapter Three. Optimization will be discussed in Chapter 
Four along with a discussion of the Automated Design 
Synthesis (ADS) program used as the optimizer routine for 
the pc)e placement and robustness design procedure developed 
in this thesis. The thesis methodology will be discussed and 
outlined in Chapter Five. This will be followed by chapters 
discussing results from selected problems. Conclusions will 
be presented in the final chapter. 
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II. SINGLE-INPUT SI NGLE-OOTPtJ? SYSTEMS 

Ihe purpose of this thesis research has been to develop 
a method of obtaining a robust multivariable control system 
design. A brief review of the concept of robustness and 
stability in the framework of a conventional SISO system 
will be dene before pursuing the concepts in a more 
complicated fashion in the following chapters. A simple 
interpretation of robustness is the ability of the system to 
tolerate design perturbations. These perturbations could be 
in the form of actuator failures, plant parameter uncer- 
tainty, unmodeled dynamics or nonlinear terms, or any one of 
many other perturbations to the nominal design of the 
system. 

The primary reason for feedback systems is the control 
cf uncertainty within the system. By appropriate use of 
feedback, properties that would lead to an unstable system 
may be controlled. When stability and robustness aspects 
are considered for a SISO system, frequency domain design 
concepts, using either Nyguist or Bode plots, are normally 
used. Robustness in SISO systems is formulated naturally by 
the concept of gain and phase margins, both of which are 
readily available on the Nyguist or 3ode diagram. 

In figure 2.1a nominal feedback system can be seen with 
a perturbation element cx! (s) placed in series with the 
nominal system. When =1 the system is nominal and stable. 
To determine the positive phase margin of the system the 
value of o(. (s) = o( { jw) = e 7 will be changed oj varying*/ 
until the system just becomes unstable. This value of (/ will 
then be the system phase margin. The negative phase margin 
can be computed in the same way. To find the gain margin the 
magnitude value of is increased until the system just 
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Figure 2.1 Basic Gain and Phase Margin Determination Model. 

becomes unstable. This is the upward gain margin. A similar 
lower margin is also defined. 

These gain and phase margins may not be adequate meas- 
ures cf robustness [ Eef . 2] because they do not account for 
simultaneous variation in both gain and phase. Therefore, 
while large individual gain or phase changes may not desta- 
bilize the system, small simultaneous changes in gain and 
phase may destabilize the system. This is not a major diffi- 
culty in classical SISO techniques because the effect can be 
easily detected. 

Gain and phase margin can be defined in terms of the 
cpen-lcop frequency domain plots in either the Bode or 
Nyguist format. Figure 2. 2 depicts a classical Bode plot 
showing gain and phase margin determination from the plot. 
The Nyguist plot may also be used to obtain this informa- 
tion. Nyguist criterion states that if the open-loop 
transfer function G (s) H (s) has n poles in the right half 
plane and the limit cf G (s) H (s) =constan t as s->*® then for a 
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Figure 2.2 Classical Bode Plot. 

stable system the locus of G(s)H{s) will encircle the -1+jO 
point n times in the counterclockwise direction as s varies 
along the Nyguist ccntour. If there are no poles in the 
right half s plane then the locus will not encircle the -1 * 
jO point. The diagram in figure 2.3 illustrates a nominally 
stable system. The gain and phase margin may be determined 
directly from the diagram. 

Any change in the loop transfer function, provided the 
order of G (s)H (s) does not change, that changes the number 
of times the locus of G(s)H(s) encircles the (-1,0) point in 
the Nyguist plot causes the system to become unstable. This 
leads to the conclusion that the minimum distance of the 
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Figure 2.3 Hyquist Plot of Stable System. 

locus of G(s)H(s) to the (-1/0) point is a measure of the 
system stability. This distance concept carries over 
directly to the MlttC system as will be shown in the next 
chapter. Examples cf a multiplicative perturbation and an 
additive perturbation illustrate this idea. Figure 2.4 is an 
additively perturbed system. Figure 2.5 shows the Nyguist 
plot for this system. Assuming that the plant is itself 
stable and the perturbations are also stable the diagram may 
then be used to determine how near the system is to insta- 
bility for the given perturbation. 
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Figure 2.4 Additively Perturbed System. 

Since the system is stable the (-1,0) point is encircled 
the correct number of times by the nominal plant. If the 
locus of g (job ) in the diagram is warped until it passes 
beyond the (-1,0) point then clearly the number of encircle- 
ments of this point will change and the system will become 
unstable, assuming the order of the plant is not changed by 
the perturbation. To keep the locus of points from moving 
beyond the (-1,0) point equation 2.1 must hold. 
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I A9 (jw) \<\ 1 ♦ g (joJ) | (2.1) 

This condition is illustrated in figure 2.6. The right-hand 
side of equation 2.1 is just the magnitude of the return 
difference transfer function of the nominal system. The 
multiplicative case is depicted in figure 2.7 with its asso- 
ciated Nyquist plot in figure 2.8. The requirement for 
stability is similar to the additive case and may be stated 
in equation 2.2 

1 &9 (jw) I <1 1 + (g (jou)-M | (2.2) 
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Figure 2.6 Hyquist for Inequality Additive Condition. 

Ihe atove arguments will be applied again in Chapter 3 to 
develop multivariable stability and robustness properties. 

Ihe linear quadratic design has been the primary method 
employed in modern control design practice. In this method 
an optimal state feedback control law is developed to find a 
set of feedback gains that optimizes a chosen performance 
index. Ihe performance index for the steady-state case is 
given in equation 2.3 
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Figure 2.7 Multiplicative System. 



PI = J(x , Qx+u’Ru)dt (2.3) 

where the x’x term and the u‘u term form quadratics. The Q 
and E matrices are chosen by the designer to provide the 
test compromise between the minimum error of the system and 
the minimum energy needed to control the system. The LQ 
method is based on the use of closed-loop state variable 
feedback fcr the control of the system. In the MIMO problem 
IQ methods have been used extensively because of their guar- 
anteed stability margins with diagonal weighting matrices. 
For diagonal weighting matrices the LQ method yields a 
guaranteed phase margin of 60 degrees and -6 db to infinite 
gain margin. 
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Figure 2.8 Nyguist Plot for Multiplicative System. 

For SISO linear quadratic optimal regulators these 
stability margins can be developed from inequality 2.4 
[Bef . 3] 

| 1 4 f * (j<JI - a) “itl > 1 (2.4) 

\ 



22 



Writing the Kalman eguation as 2.5 



|1 4 f « (sl-a)-ib| 2 = 1 + (1/^)1 G ^ { s) | 2 (2.5) 

cne has further 

| 1 + G (S)H(S) | 2 = 1 + (1/f) IG K (S) | 2 (2.6) 

and for all s = j c*j and 0< co < o° the function (1 / p) G* (s) is 

greater than zero, therefore the Kalman inequality is shown 
to be 2.7. 

| 1 4 G(s)H(s) | > 1 (2.7) 
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The graphical interpretation of this result is that the 
polar plot cf G(s)H(s) must remain outside the unit circle 
centered at the -1 + jO point for all frequencies. Figure 




Figure 2. 10 Polar Plot of an Optimal System. 

2.10 shows such a polar plot. Since the optimal regulator 
with the loop broken at the input to the plant (denoted by 
the x in figure 2.9 ) does not penetrate the unit disk about 
-1 + jO this means that the single input regulator will have 
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a phase margin of at least 60 degrees and a gain margin 
tolerance of fifty percent gain reduction and infinite 
upward margin. Further discussion of this property may be 
found in Anderson and Iloore [Ref. 4], With this hasic 
review of the concepts of stability and robustness in the 
classical SISO system complete, the next chapter will extend 
some cf these basic concepts to the HIMO system. 
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III. MOLTI VARIABLE SYSTEMS 



linear quadratic design has performed relatively well in 
aircraft control concepts because of the ability tc formu- 
late the system state equation and to quantify acceptable 
performance indices for the system. Industrial applications 
for L£ theory have been less successful and have led British 
researchers to look into forms of decoupled design methods. 
Cne of the primary methods used in multivariable design is 
to make the system totally decoupled. This method allows 
each loop tc be designed as a separate entity by classical 
means. Cne of the primary difficulties with this method is 
the problem of finding a compensator which will totally 
decouple the system. The method also suffers from the 
effects cf cross-coupled perturbation terms. A method that 
does' not totally decouple the system but makes the design 
problem simplier by designing a compensator that only causes 
the diagonal terms of the transfer matrix to be dominant 
over all cff-diagonal terms has also been developed, 
[Sef. 5]. Classical frequency domain techniques are then 
used to design each loop of the system. The major difficulty 
with this and other single loop design techniques is their 
failure to account for cross-coupling perturbation terms 
that may interact between the loops. The individual loops 
may be highly robust in these designs but the overall system 
robustness may be lew because of the loop interaction not 
accounted for in the design. This is the precise area that 
the singular value analysis has proven so beneficial in 
design . Singular value concepts may be applied tc conven- 
tionally designed systems to assess their robustness. Bor 
instance, if a system is designed by LQ methods the designer 
may then formulate the transfer function of the system and 
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assess tie singular values ox the return difference matrix 
of the system. If the minimum singular value is found to be 
low at some critical frequencies the designee can then 
modify the Q and/or E matrices chosen in the LQ performance 
index and recalculate the design. In this iterative fashion 
a rotust design would be developed. 

A generalization of the SISO Nyguist theory discussed in 
the previous chapter has been made for the MIHO problem. 
This generalization leads directly to the singular value 
concept. The generalization is expressed in the form of the 
multivariable Nyguist theorem which requires that a closed 
loop stable system have the same number of counterclockwise 
encirclements of the origin by the locus of the det (I+G (jen) } 
as the number of open loop poles that are unstable. This 
theorem is formally stated as; 

let N[f (s) ] denote the number of clockwise encirclements of 
(-1,0) by the locus of f(s) as s traverses the contour D 
of figure 3.1 in a clockwise sense. The closed-loop system 
will he stable if and only if for all R sufficiently large 
N[f(s)] = -P 

where 

f (s) = -1 + det[ I + G (s) ] = </V(s)/ ^ cl (s ) -1 and 

P = the number of closed right-half plaa e zeros of 
(£) . 

The application of the Nyguist theorem comes through the 
fact that a multivariable system will not be robust to 
modelling errors if the return difference matrix, I + §, is 
nearly singular for some frequency. If I + G is nearly 
singular a small change in G may make I + G exactly 
singular. This causes the det (I + G) to become zero and the 
Nyguist encirclement count to change indicating an unstable 
system. It is possible for very small changes in I+G to 
produce large changes in the determinant of I+G. The matrix 
I+G 
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Figure 3.1 Nyguist D Contour. 

«■ 

[ 10/s+a 9. 99/s+a 

10/s+a 10/s+a 

has determinant 0. 1/ (s+a) 2 . If the element p i( is changed by 
only cne percent to 9. 9/s+a the determinant becomes 
1.1/ (s+a) 2 which is a significant change in the determinant 
value. Therefore, it is evident that det (I + G| is net an 
accurate measure of how near the return difference is to 
singularity. Researchers, in the field of controls 
[Ref. 6], [Ref. 7 ], [Ref. 8], [Ref. 9] have used singular 
value analysis to determine how near the return difference 
matrix is to singularity. 

Since the number cf encirclements of the Nyguist diagram 
changes as f (s) passes through the -1 point or when det (I+G) 
is zero it is important to find how near the return differ- 
ence matrix I+G is to being singular. This nearness to 
singularity can be interpreted as closeness of the matrix G 
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to the critical point, -1. A quantity which can be used to 
express the nearness to singularity of the matrix is the 
minimum matrix singular value denoted by . Given a matrix 

A the singular value may be expressed by equation 3. 1 

C (A) = min ( U U A)) /Z (3.1) 

where \ (£*A) is the eigenvalue of the complex conjugate 
transpose of A times A. A basic MI MO linear system is 




Figure 3-2 Basic Multi-input, Multi-output System. 

depicted in figure 3.2. An additive perturbation to the 
plant is shown in figure 3. 3. If the plant is stable before 
the perturbation is added to the system the Nyquist theorem 
will be satisfied and the locus of GH will not encircle the 
-1,0 critical point. When the perturbation is added to the 
system as long as the Nyguist locus is not forced to 
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Figure 3.3 additive Perturbation. 

encircle the -1 point the system will remain stable. A 
sufficient condition, recalling the S ISO discussion in 
chapter 2, for the perturbed Nyquist plot not to change 
encirclements is that the norm of the perturbation &G 
remain less than the norm of the return difference matrix as 
expressed in equation 3.2. 

I I ^ £(j*») ||< 1/1 ! (1+5) -») | j (3.2) 

>0 This condition will guarantee that the Locus of the 
det (I+G) does not pass through the -1 point. If the 1^ or 
Euclidean norm is assumed for this condition the equation 
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3.2 may be expressed in terms of singular values as equation 

3.3. 

ff(AG)< ^(1 + G) (3.3) 

This result states that as long as the maximum singular 
value of the perturbation matrix AG is below the minimum 
norm value cf the return difference matrix the system will 
remain stable. The problem of guaranteeing robustness 
becomes that of finding the largest norm of the perturbation 
quantity, the largest singular value, for which the smallest 
norm or singular value of the return difference matrix will 
remain non singular. 

The multiplicative form for a system such as figure 3.4 
gives the similar norm equation in 3.4 

I I Aq,(jw) | J < 1/1 1 (i + .(G)- M" 1 1 1 (3.4) 

OJ >0. which may be expressed as; 

CT { AG) < <d ( I ♦ G-i ) (3.5) 

Singular value decomposition software is readily available 
to determine how near the matrix I+G or 1+ (G)“ l is to 

N. /V -V 

singularity. 
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Figure 3.4 Hulti plica tive Perturbation. 
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IV. OPTIMIZATION 



The purpose of this chapter will be to briefly describe 
several of the currently employed optimization techniques 
and the Automated Design Synthesis (ADS) program, [Ref. 10], 
which employs these techniques. In general, optimization 
implies finding the "best" possible solution to a problem. 
In actuality the best solution found by an optimization 
technigue might really only be a "better" solution to the 
problem. The purpose of ADS and other optimization routines 
is to allow a rational search to be conducted to find the 
best possible design. The techniques of numerical optimiza- 
tion are used to logically vary the various parameters that 
affect the design until a good solution is found. 

As an example of an unconstrained optimization problem 
consider the following problem developed in [Ref. 11]- The 
problem is to minimize the function 

F (x) = 10x y *-20 x, 2x z + 10x £ 2+x / 2 -2x, +5 (4.1) 

F (x) is often called the objective function, the cost func- 
tion or the penalty function. Since there are no conditions 
imposed on the design variables, x, and x and no addi- 

tional limits imposed on the overall design, the problem is 
considered to be one of unconstrained minimization. Figure 
4.1 represents this problem in the design space. From the 
figure it appears the optimum is near the point 1,1. 
Calculus may be applied to determine the optimum exactly. 
Taking the derivatives 
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Figure 4. 1 Design Space for Exaaple Problem. 



pF ( 2 )/^x, =40x ( 3-40X, x z +2x, -2=0 (4.2) 

9 ? (x) / 9 x t =-20x, 2 + 2 0x^=0 (4.3) 

and then solving the set of eguations it is found that x, 
=1 . 0 and x x =1 . 0 . 

If design conditions are imposed on the problem then the 
optimi2ation becomes one of constrained function minimiza- 
tion. In other words, while the minimum of a function is 
still sought, this minimum must exist within the limits 
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imposed by the design conditions or constraints. The minimum 
of an unconstrained function will not necessarily be the 
same minimum for a constrained function. 

As another example , also from [Ref. 11], a design is 
sought which gives the minimum weight of a particular 
column. This weight is expressed as: 

B= y>Ah= jo/Z Dh (4.4) 

where is the unit weight of the material and A the cross 
sectional area. The stress in the column is given by: 

CT =P/A = ?/* Dt (4.5) 

without going into detail the design will be constrained by 
the allowable stress on the structure. Other constraints , 
Euler buckling and shell buckling, are also of interest to 
the designer of this column. The design problem is then 
stated as: 

minimize W = Eth (4.6) 

for the constraints 



g (1) =cr /c r - 1< 0 



(4.7) 



g (2) = <3 * /(To ~ 1 ^ 0 



(4.8) 



9 (3) = Gy*- -1< 0 



(4.9) 



with B> 10e-06 and t >10e-0 6. 
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g (4) = t-D< 0 



(4. 10) 






& *C fj ' &s refer to the allowable design stress, Euler 

buckling stress and the shell buckling stress respectively.. 




Figure 4.2 Design Space for Column Problem. 

Figure 4.2 is a design space diagram for this problem. It 
is noted from the figure 4. 2 that the optimum is not unigue 
and can be any value along the arc noted as optimum in the 
figure. 
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The general optimization problem is then written as: 
minimize F(x) objective (4.11) 



subject to g • (x) < 0 j=1 y m inequality constraint (4.12) 
0 



h K (x) = 0 k=1 ,1 equality constraint (4.13) 



x^< x.< x'l' i=1,n side constraint (4.14) 

( l L 

where X = col (x1,x2,.. . ,xn) is the design variaDle. 

The methods used to solve this problem are usually 
iterative. After the establishment of an initial set of 
variables the optimizer will update this initial value until 
the optimum values are found. Again /borrowing from 
[Ref. 11], the iterative technique may be demonstrated by a 
simple example. Figure 4.3 is used to illustrate this 
problem. Given the initial data set X° the formula 

XI = X« + <x? S 1 (4. 15) 

can be used to upgrade this estimate of X. The vector S is 
the search direction for the iteration and the scalar quan- 
tity X is the distance of the move in the S direction. 
Beginning at X° it is desired to reduce the objective func- 
tion. The search for values of X that reduce the objective 
function is made in the $ direction which in this example is 
the opposite of the gradient of the function at point X°. 
The choice of S could be arbitrary as long as it reduces the 
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Figure 4.3 Illustrative Example for Iteration. 

value of the objective. If S is chosen to be opposite the 
gradient of the objective then the search would be a 
steepest descent search. Once the direction is chosen then 
the scalar must be found to minimize the objective along 
this direction vector without crossing a constraint 
boundary. The technique for finding e< is to evaluate X, the 
objective and the constraint functions for values of using 
computer code and numerical interpolation to estimate <x . 
This value then gives the minimum value for F (X ) along 

this search vector. Now that c< for this direction is known 
the new equation becomes: 

X 2 = x 1 ♦ **S2 {4. 16) 

and now X 1 is used to start a new search in direction § 2 
and compute a new <* # that continues to reduce the objective 
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function. These one dimensional searches are continued until 
no more reduction in the objective can be found. At this 
point no further design improvement is possible. 

Optimization techniques do not always lead to the abso- 
lute optimum when applied to problems of practical interest. 
The reasons for this could be numerical ill-conditioning of 
the problem formulation or simply that there are multiple 
solutions to the problem. Because of these difficulties it 
may be adviseable to choose several different starting 
points for the optimization process and use engineering 
judgement as to the design most applicable to the problem 
under analysis. 

Considering the unconstrained case first where the 
desire is to minimize the function F (x) , it is well known 
that F (x) will have a minimum where the gradient of F (x) is 
zero. That is: 

grad (F (x) ) =0 (4.17) 

with the y F (x) defined as: 

grad (F(x)) = ( 2F (xJ/Qx, , . . . , 3F(x) /3x„) ' (4.13) 

Figure 4.4 shows why this is a necessary condition but dees 
not guarantee a global minimum. The gradient of F (x) is 
zero at all four points A,B,C,and D. However, only A and D 
are minima. A would be the global minimum for the function 
as defined here. D would be only a relative minimum. To 
check that the zero gradient corresponds to a minimum the 
Hessian matrix, i.e the matrix of second partial derivatives, 
can be examined for positive definiteness. A positive defi- 
nite Hessian ensures a relative minimum. The only way to 
prove a global minimum for the function is to show that the 
Hessian matrix is positive definite for all design variables 
X. A test that is seldom possible to perform. 
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If the problem of minimization is a constrained problem 
the situation is different from that discussed above. The 
objective function gradient does not have to be zero at the 
optimum. Figure 4.5 illustrates this case. Using figure 4.5 
and assuming a start at point A it is necessary to choose a 
search direction S which will reduce the objective function 
while not violating the constraint functions. Any direction 
that will reduce the objective is said to be useable. This 
is seen to be the half plane sector below the tangent to 
F (z) at point A. If a half plane to the right of the 
tangent to the active constraint at point A is considered 
then the feasible sector is formed. A combination of the two 
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conditions gives the useable, feasible sector where the 
search direction S most be chosen. In mathematical nomencla- 
ture the above argument may be stated as; 

useable direction grad { F (X) ) • S< 0 (4.19) 



feasible direction grad ( g_. (x) ) • S< 0 for all j (4.20) 
for which g ^ (x) = 0 . 

Feint 3 in figure 4.5 shows a point where the gradient 
of the objective and constraint point in exactly the oppo- 
site direction. At a point such as 3 the only search vector 
S that meets requirements for useability and feasibility is 
tangent to the constraint boundary and to a line of constant 
objective function. This condition is stated as: 

grad (F (x) ) +Z>igrad (gj (x) ) + 2>^ t grad(h x (x) ) = 0 (4.21) 

where > 0 and ^ -unrestricted. 

Kith this brief background in optimization complete the 
Automated Design Synthesis (ADS) program will be briefly 
discussed. More detail on the ADS routine can be found in 
[Ref. 10 ]. This code was developed as a follow-on to the 
successful CONMIN code [Ref. 12] developed by Va nderplaats. 
It is designed as a black box optimizer which allows the 
user to choose combinations of one dimensional search, opti- 
mization algorithm and optimization strategy. These will be 
discussed later. For the user with specific requirements the 
code may be tailored by parameter modification to meet 
specific requirements. For the work done in this thesis it 
was assumed the user of the code has no detailed knowledge 
of optimization and will want to use the code in the 
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Figure 4.5 Constrained Optimization Example. 

simplest form. As such, the ADS calls made from the main 
program use only default parameters and first focward finite 
difference gradients. Should analytical gradients be avail- 
able the user could use them within the code if desired with 
no difficulty. The calls made by the user to the ADS 
routines specify several important aspects of the problem 
solution. The user may call for an optimization strategy to 
be used in the routine. This is not required and its use 
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depends cn the problem. These strategies are discussed in 
£Bef. 10]. Two strategies used most often in the analysis 
done for this thesis have been seguential unconstrained 
minimization using quadratic exterior penalty function and 
the augmented Lagrange multiplier (AIM) method. Others 
available include seguential linear programming and seguen- 
tial guadratic programming. 

The basic optimizer is also chosen by the user from two 
unconstrained and three constrained optimization algorithms. 
The unconstrained algorithms are Fletcher-Hee ves conjugate 
directions, Davidon-Fletcher-Powell (DFP) variable metric 
method and the Broydcr.-Flet cher-Goldf arb-Shanno (BFGS) vari- 
able metric method. The method of feasible directions and 
robust feasible directions are available for constrained 
minimization. 

The user has available several types of one dimensional 
searches using Golden Section or polynomial approximation 
technigues. The ADS code has tailored these one dimensional 
search algorithms for the unconstrained and constrained 
cases, allowing the user to make appropriate choices for the 
type of problem to be solved. Figure 4.6 shows the basic 
organization of the AES program. 

Since it is not the purpose of this chapter to cover the 
subject of optimization in detail all the possible routines 
in ADS will not be discussed. A few of the routines found to 
work well for the work completed in this thesis will be 
briefly described, however. There are several methods of 
optimizing functions of one variable or one dimensional 
searches. For instance a large number of points, n, could be 
chosen and the function F(x) evaluated at each point. The 
point corresponding to the smallest value of F (x) could then 
be considered the optimum value of the one dimensional 
search. This method is hit and miss and better methods of 
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USER-SUPPLIED CALLING PROGRAM 




Figure 4.6 Organization of ADS Prograa. 

locating optimum points can be used. The methods used in ADS 
are Golden Section and polynomial interpolation. 

Golden Section search methods are easy to program on the 
digital computer and do not require continuous derivatives. 
They tave a known convergence rate and are reliable for 
poorly conditioned problems. The major drawback of the 
Golden Section routines is the large number of function 
evaluations require. The Golden Section algorithm is simply 
illustrated through the use of figure 4.7 Assume that X° and 
X 1 are known to be bounds on the curve’s minimum value. 
Also, the function values of F (x) , F° and F l ,are evaluated 
and known at these points. By picking two intermediate 
points X 2 and X 3 where X 2 < X 3 and evaluating the function 
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Figure 4.7 Golden Section Diagram. 

at these points a hound on the minimum may be modified. 
Since F 2 for this figure is larger than F 3 the point X 2 
forms the new lower bound. The minimum is now between X 2 and 
X 1 . If the function F 4 at X 4 is determined and shown to be 
larger than F 3 then X 4 becomes the new upper bound on the 
minimum. Ey repeating this procedure the bounds may be 
narrowed to any desired tolerance. The Golden Section rule 
is applied to this problem to reduce the bounds in the 
quickest possible time. By appropriately picking the values 
of X 1 s at which each function evaluation is made an effi- 
cient algorithm that uses the ratio proportion of the Gclden 
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Section, i.e. X 2 /X l = 1.62803 is developed. By choosing a 

value 'T based on this Golden Section rule where ^ = 0.38 197 
estimates for interior points X 2 and X 3 can be made as: 

X 2 = (1 - £ ) X° + 2-Xi (4. 22) 



X 3 = 't X° + (1 - 't) Xi (4. 23) 

As each new bound is found the process repeats until the 
accepted level of convergence is reached. 

The polynomial interpolation method is accomplished by 
first fittirg a polynomial curve to the points about where 
the minimum is desired and then finding the minimum of the 
polynomial function. For example if the function F(X) is 
approximated by a quadratic as: 

F = a 0 +a ( X + a z X 2 (4. 24) 

sc 

Then the value of X, X , where F* is zero can be shewn 
to be : 

X* = -a , /2a z (4. 25) 

If a ^ is positive then F will be minimum. Other degrees of 
polynomials may be used in similar fashion. 

Now that the basic one dimensional search methods have 
been reviewed the next step is to examine the basic opti- 
mization routines. First, the unconstrained case will be 
reviewed. The optimum X*is at the point where 
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grad ( F </) ) = 0 



(4. 26) 



Several zero order methods exist for the purpose of finding 
the minimum value. These include random search, Powell’s 
method and’ Box’s method. Since ADS does not use these 
methods they will not he discussed. 

Automated Design Synthesis makes use of first order 
methods which will new be discussed. The steepest descent 
method is best known, but poor in performance. Just as the 
name implies a search direction is chosen opposite the 
gradient of the objective function. At iteration 0 

SO = - g r ad ( F (X°) ) (4. 27) 

Figure 4.8 shows this algorithm geometrically. Note that 
the method simply stair steps its way down the "hill” to the 
valley or minimum. AES uses the Fletcher-Reeves modification 
to steepest descent. In this routine a conjugate direction 
is chosen to improve the speed of the search. Figure 4.9 

shows how this method tracks to the solution. 

The variable metric methods listed earlier are usually 
more powerful than Fletcher-Reeves because they store infor- 
mation that allows tie algorithm to approximate the inverse 
of the Hessian matrix or second derivative. For further 
discussion of these methods see [Ref. 11] or other similar 
optimization texts. 

AES employs two direct methods for constrains d minimiza- 
tion. Cne method is that of feasible direction and the ether 
is the method of robust feasible direction. Since these two 
methods were seldom employed within the work presented in 
this paper they will not be discussed. 

The methods chosen to handle the constrained minimiza- 
tion problems formulated in this thesis are referred to as 
Sequential Unconstrained Minimization Techniques (SUMT) . 
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Figure 4.8 Steepest Descent Algorithn. 

SOMT methods are methods which formulate the objective 
function and the constraint functions into an augmented 
objective function and then solving the problem as if it 
were an unconstrained optimization task. ADS employs 
interior and exterior penalty function techniques as well as 
an Augmented Lagrange Multiplier (AIM) technique. 

The exterior penalty function method is incorporated by 
forming a penalty from the constraint equations. This 
penalty is of the form : 

P(x) =21 max (0 # g. (x))2 + -Sl (h (x)) 2 (4.28) 

J*i i K-*> 



48 



*2 




Figure 4.9 Fletcher-Reeves Conjugate Direction Algoritha. 

From this eguation it is obvious that P (x) is zero if all 
the constraints are satisfied. That is, if g . (x) < 0 and \ 

(x) = 0, then all conditions of the penalty function are 
satisfied. If an element of the penalty function is violated 
then the penalty increases as the square of the violated 
constraint. A pseudo or augmented objective function is 
formulated where : 

F (¥) + P (x) (4.29) 

The constant r^ is a weighting factor for the penalty. It is 
adjusted with ADS as the optimization proceeds to allow the 
program to systematically converge to an optimum solution. 
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Cne disadvantage of this method is that if the optimizer is 
stopped short of the optimum the design will be in the 
infeasible region and probably not useful. 

Originally developed as a method to solve equality 
constrained problems the version of AIM in ADS has been 
modified to work with both equality and inequality 

constraints. The statement is as follows: 

minimize ? (x) (4.30) 



subject to h K (x) = 0 k=1 ,1 

Next a lagrangian is created such that 



I (x,» 



F (X) 



zK 



h*(X) 



(4.31) 



(4.32) 



Then the equation is augmented with an exterior penalty 
function such that 



A{x /( k,r p ) = F (X) +£(> lt h K (x) + r p (h K .(x))2 ) (4.33) 

The power of this method is that in theory precise 
constraint matching is possible whereas in the exterior 
penalty function methcd it is not. The full details of the 
method will not be covered here, however, the final form of 
the objective function will be included: 




) (4. 34) 



where =max (g j (x) , 2r P ) , and the update formulas for 

the Lagrange multipliers are: 
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( 4 . 35 ) 



/j + '=A^ + 2r r (max (g^ (x) # ->j/2r p ) ) j=1,m 

C. = >L* 2t p h i<*> C- 36 ) 

Kith this trief review of optimization and the ADS program 
completed the next chapter will discuss the program develop- 
ment for this thesis research- 
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V. OPTIMIZATION DESIGN PROCEDURE 

Tc accomplish the task of designing a control system 
with acceptable time domain performance and robustness 
characteristics in a straight forward manner a numerical 
optimization technique was chosen as the method of iiplemen- 
tation of the design algorithm. Using numerical procedures 
to adjust selected design variables, in this case the feed- 
back and/or filter gains, a desired level of performance can 
he achieved. This level of performance is actually a combi- 
nation of time domain performance and robustness or 
frequency domain performance. 3y establishing the criteria 
for the system performance in terms of an optimization 
objective and constraint functions a versatile procedure can 
he developed to set the system feedback gains and affect an 
acceptable design in terms of performance and robustness. 
The pcle placement and robustness (POPLAR) design procedure 
uses pole placement to establish a designer selected 
performance level and then a minimum singular value level to 
establish the robustness. 

The pole placement portion of the procedure will be 
discussed first. The pole placement technique was chosen 
because it was relatively easy to implement through a numer- 
ical optimization routine. By using this numerical procedure 
it is also simple to incorporate robustness into the proce- 
dure along with the performance requirements. A numerical 
technique similar to one posed in £Ref. 13] was chosen for 
the pcle placement algorithm. An unconstrained optimization 
routine from the IBM IMSL library was used for this program. 
The routine, a Newton method, uses adjustments to the output 
feedback gains to reduce the size of an objective function. 
This objective function was expressed as a function of the 
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pole location of the system and as the objective was reduced 
the poles were moved toward the desired locations. To 
provide more versatility in the pole placement algorithm a 
method that can use constraint functions as well as uncon- 
strained optimization was chosen for this program. The 
designer may use either objective, constraint or a combina- 
tion of functions to secure the desired pole locations. As 
currently implemented in the program the cost or objective 
portion of the pole placement procedure is constructed as 
equation 5. 1 

OBJ = r (A*,.- W; ) 2 + (Ax*. -Ai; ) 2 (5.1) 

where Ar. = real eigenvalue 

)*i = imaginary eigenvalue 

= desired eigenvalue location 
= desired eigenvalue location 
The constraint formulation is a function tiat must be 
kept negative or the constraint is violated. It is written 
as equation 5.2 

g(j) = V ) 2 + )2 - r (5.2) 

where r is a tolerance circle established as a function of 
pole placement position. Since the aim of the optimizer is 
to keep g negative any time the A function of the 
constraint is greater than r the constraint will become 
active, i. e. violated. The optimizer will then attempt to 
move the constraint to the inactive status by adjusting the 
design parameters of the system. 

Consideration of implementation of the frequency domain 
or robustness portion of the design procedure begins with 
the concept of MIMO phase and gain margins. Several useful 
theorems on singular value analysis of multiloop systems are 
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presented in [Ref. 8]. One of these theorems relates the 
matrix singular value of the return difference function to a 
parameter, o( , and further shows that as long as the maximum 
singular value of the perturbation function (L - 1 - I) 

remains less than this o(, the system remains stable. The 
value of o( is then related to gain and phase margins of the 
MIMO system. The relationship developed is given in equa- 
tions 5.3 and 5. 4: 

gain margin = GM = 1/(1+£><o) (5.3) 



phase margin = PM = ±cos " l ( 1-<* 0 2 /2) (5.4) 

provided that equation 5.5 holds. 

pMI + G) > Ko 
for seme <><„< 1 

These phase and gain margins are guaranteed in 
simultaneou sly . 

Universal gain and phase margin curves, [ Hef 
on the minimum singular values of the return 
matrix are developed from equation 5.6. 

(I-i-I) =aaxy (1- 1/k A ) 2 + 2/k rt ( 1-cos (5.6) 

for all n with k„ > 0. These curves shown in figure 5.1 
allow the designer tc pick a singular value that corresponds 
to a specific gain and phase margin for a given system. In 
addition tc the universal gain and phase plot [Bef. 15] 
developes an optimizer solution for formulating a robust 
controller using the CONMIN algorithm [Ref. 12]. 



(5.5) 

every loop 

. 14], based 

difference 
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Figure 5-1 Universal Gain and Phase Singular falue Plot. 

Since the universal curve in figure 5.1 provides a 
convenient method of specifying gain and phase margins in 
terms of singular values the robustness portion of the pole 
placement and robustness design procedure uses the minimum 
singular value level of the return difference matrix to 
determine the robustness. The minimum singular value level 
is formulated as an objective or constraint function in 
equation 5. 7 

J=Z.(max (0, (<£-«- (j»,p)) ) 2 (5.7) 

J 

The optimization procedure may be used to change feedback 
gains until the minimum singular value is raised above this 
desired design level. Although the same formulation can be 
used as a negative constraint function it has not been 
implemented as such within this program. There are numerous 
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ways the singular value formulation could be implemented 
within the program by changes of the code if design require- 
ments forced such changes. 

The pcle placement and robustness design program is 
based on the ADS code discussed in Chapter 4 to implement 
the design variable selection procedures. The pole placement 
and robustness program consists of two separate programs. 
The first program is used to provide designs for state or 
output feedback problems while the second program is used 
for observer or filter designs. 

For the state or output feedback design program the user 
must input the plant matrices A,3,C and initial starting 
values for the feedback matrix F. The matrices correspond to 
the following linear differential system: 

x = Ax + Bu (5.8) 



y = Cx 

Z- — 



(5.9) 



u = -Fx (5. 10) 

A feed-forward matrix has not been considered in the current 
program. A feed forward matrix could be added to the proce- 
dure if required for specific design cases. 

As the design program is currently coded the user may 
run output feedback or state feedback design by specifying 
the C matrix as the diagonal (I) matrix for state feedback. 
The program relies on initial starting values of the feed- 
back gains, F. As discussed in Chapter 4 there is no 
guarantee that the optimum found by the procedure each time 
is the global optimum or that the procedure will always 
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converge to an acceptable solution. The ability to select 
acceptable starting values for the feedback gains will make 
the procedure more efficient in operation. As currently 
employed, the program is used to obtain pole placement and 
robustness for a given set of starting gains and a selected 
optimization routine from the ADS program. If the optimizer 
is net able to meet the desired design goals on this program 
run two options are available. First, change to a different 
optimization routine from the list of available ADS routines 
and rerun the problem. This was usually successful in 
improving the design. Second, the designer uses a new set of 
starting values for the feedback gains and repeats the 
design procedure. Both options might be used on particularly 
difficult cases. 

The pole placement and robustness design procedure has 
consistently been able to find improved designs; however the 
program does not always yield acceptable designs. Certain 
problems reguire changes in the optimizer routine and modi- 
fication in the initial feedback gain starting values in 
order to obtain acceptable designs. Using the IBM 3033 time 
share system the pole placement and robustness routine 
reguires about ten CPU seconds to work a second order 
problem and on the order of 15 to 60 CPU seconds to run a 
forth order problem. The actual amount of time varies with 
optimization reguirements and time share utilization. 

The observer robustness design program reguires two 
passes of the ADS program. In the first pass the feedback 
gains, F, of the controller are computed to obtain the 
desired pole locations. The second pass of the ADS routine 
is used to adjust the observer gains to recover the system 
robustness. The two pass procedure was chosen because it 
allows a smaller number of design variables at each stage of 
the optimization and much more efficient computer usage. 
Figure 5.2 shows how the observer is implemented. This 
diagram is algebraicly stated as; 
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where x is the state, x is the estimator variable, F is the 
feedback gain and K is the observer gain. The design cf the 
feedback gains and the observer gains are accomplished as 
separate quantities in keeping with the separation Frin- 




Figure 5.2 Observer Implementation. 

ciple. In using tie pole placement and robustness design 
procedure for the observer system, initial values of the F 
and K matrices must be input. The same or different opti- 
mization techniques from ADS may be employed. 

The routines contained in this thesis have been based on 
the input additive singular value level. The pole placement 
and robustness design algorithm computes input additive. 
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output additive, input multiplicative, and output multipli- 
cative singular values. Any of these singular values can be 
incorporated into the objective or constraint formulations 
but have net been for this version of tne program. 

In summary, the pole placement and robustness design 
procedure is a straight forward numerical optimization 
procedure fer the practical application modern MIMO system 
analysis. The new aspects of the procedure are the implemen- 
tation of both pole placement and robustness criteria within 
the same design program. The versatility of the pole place- 
ment and robustness design is obtained by incorporating a 
state of the art optimizer routine ADS, with currently 
available singular value computation routines. Using the 
optimizer format for the pole placement and robustness 
design gives the designer the ability to modify variables 
directly that affect both the time domain or performance of 
the system and the frequency domain or robustness of the 
system. The numerical optimization incorporated into the 
pole placement and robustness design is flexible enough to 
incorporate the double pass design technique using the sepa- 
ration of the feedback and filter gains to obtain robustness 
recovery for observer based controllers. 
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VI. INTRO EOCTORY PROBLEM 



The purpose of this introductory problem is to review 
some results of modern multivariable robustness theory using 
a simple problem. This same simple system is then used as a 
test problem for the optimization technique developed for 
this thesis. The problem provides excellent insight into the 
cross-coupling problem and demonstrates how effectively the 
pole placement and icbustness design procedure can be in 
dealing with the cross perturbation terms. 

The problem chosen for this introductory analysis comes 
from £Ref. 9]. Figure 6.1 is a diagram of this basic 
system. In this problem a simple plant is specified by the 
following linear system: 
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y 1 = x 1 (6-2) 

y2 = x2 

A feedback compensation of the form of equation 6.3 was 
assumed . 
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Which gives a closed-loop system, equation 6.4 
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This system has eigenvalues at -2 
stable. Osing equation 6.5 
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and is 



(6.4) 

therefore 
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G = C(sI-S)- l £ 

the transfer matrix may be written as equation 6. 6. 

G = Fl/s+1 b , c /s + ll 

I 0 1/s+l j 



which gives the return difference matrix# equation 



I + G(s) = 



s+2/s+ 1 t /s+ 1 
0 s+2/s+ 1 



(6.5) 



( 6 . 6 ) 



6.7. 



(6.7) 




♦y 2 



Figure 6.1 Basic Multi-input Multi-output System. 



The problem 
in establishing 



shows the inadequacies of classical methods 
the robustness of the system. A brief review 
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cf the results will he presented. Using the return differ- 
ence matrix (1+5)/ the determinant may be written as 
eguation 6. 8. 

det (I+G) -1=2s + 3/ (s + 1) 2 (6.8) 

Therefore, the multivariable Nyguist diagram will be as 

Im 




shown in figure 6.2. The diagram does not encircle the 
(-1,0) point and is indicative of a closed-loop stable 
system. Considered as a SIS 0 system one has a gain margin of 
-1/3 to co and a phase margin of ±106 degrees. Under this 
criteria one can conclude that this is a good design. This 
will be shewn later net to be the case. 

Multi-input systems are often designed a loop at a time. 
Applying loop design to the system the transfer functicn for 
either lcop taken separately becomes eguation 6.9 



62 



G= 1/S + 1 



(6.9) 



figure 6.4 shows the Nyguist diagram of this problem. Using 
this Nyguist diagram the system is indicated to be stable 
and have phase and gain margins of GM=(-1,<») and 

FM=(±180). This analysis does not show the true nature of 
the robustness of the system. Since the factor b is not a 
parameter in either of the Nyguist curves it plays no part 
in the stability determination using these diagrams. 

Using the criteria of singular values discussed in 
Chapter 3 a measure of the nearness to instability for this 
problem may be obtained by plotting the minimum singular 
value of (I+G) , (T • (for numerical calculations a value of 
b /2 = 50 is assumed). Figure 6.3 shows the plot of this 
value vs. freguency. This gives a minimum singular value of 
about -23 db or 0.071 which corresponds to a gain margin of 
about 0.S3 to 1.08 and a phase margin of ±4.1 degrees. These 
phase and gain margins are guite small and are evident in 
the cross-feed perturbation problem. 

A perturbed system as shown in figure 5.5 can be 
produced which will lead to stability problems with very 
small values of perturbation. The closed-loop system if 
perturbed by a small perturbation, 5/b, L , where h lt is a 
large number, will have as a characteristic eguation, egua- 
tion 6.10 



(s I- A) = s2 + 4s +S (6.10) 

with the eigenvalues of s=-2± There is one positive 

root in this solution and the system is unstable. 

To determine the nature of the robustness of this system 
the return difference matrix must be considered. If the 
return difference matrix of the transfer function (I+G(j«)) 
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Frequency (rad/sec) 

Figure 6.3 Hiniaua Singular Value Plot, Example Problem. 

is nearly singular at some frequency co 0 , then the multivari- 
able system will not be robust with regard to any modelling 
errors within the system. This is because any small change 
in G ( joj) can then male I+jG ( jw) singular and the det (I+G) 
becomes zero, thus changing the encirclements of the Nyguist 
stability point and indicating a system instability. Using 
equation 6.1 the pole placement and robustness design method 
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Figure 6.4 Nyquist Diagram of 1/s+1. 

can be demonstrated . The first step is to establish a base- 
line for the design. The state unity feedback model in this 
simple problem gives pole locations of -2 and -2. This set 
of eigenvalues -2, -2 were chosen as the baseline for the 
system. Since a second order system reguires only two feed- 
back gains to place the poles, the diagonal feedback gains 
were chosen for pole placement purposes. The pole placement 
and robustness design program was then used to place the 
poles and in the process obtained the gains required to do 
this (1,1). A plot cf some of the singular value criteria 
obtained is shown in figure 6.6. Pole-zero plots of the 
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Figure 6.5 Perturbed System. 

closed- loop transfer functions of the closed-loop transfer 
matrix are shown in figure 6.7. In this case the poles and 
a zero are clustered about the -2 point and in the input two 
to output one channel a zero is located at the -1 point on 
the pole-zero diagram. This point closely corresponds to the 
minimum singular value frequency. The only significant 
observations are the ability of the pole placement and 
robustness routine to place the poles and the relatively 
poor singular values indicative of low robustness. 

Since there is a requirement for two feedback parameters 
to set the poles of a second order system there should be no 
additional freedom in design to account for robustness. 
Case two was a run tc demonstrate this fact. Again allowing 
only two design variables for the problem the pole placement 
and robustness desigt program was run but with an objective 
function formulated to adjust the singular value above a 
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Figure 6.6 Singular Value Plot for Simple Problem. 
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base level of 0.6. This level would correspond, using the 
universal gain margin chart, to a gain margin of -4 db to 3 
db and a phase margin of ± 35 degrees. With only the two 
design variables to work with the pole placement and robust- 
ness design program was unable to place the poles and adjust 
the singular value level to the required value. The poles 
were placed at -2, -2 but the singular value minimum was 
still cn the order of -23 db or 0.07. The pole-zerc plot 
remained almost unchanged. It is clear that additional 
degrees of freedom for the pole placement and robustness 
design program must be opened if robustness is to be 
accounted for. 

Case three was then run on the pole placement and 
robustness program by adding on additional degree of 
freedom. This case used the f gain as the additional 
design variable. A good choice as will be seen. Allowing 
the optimizer routine the extra freedom to adjust the addi- 
tional feedback gain term an excellent design was found. 
The singular value minimum became 0.885 as shown in figure 
6.8. Using the universal gain margin chart this corresponds 
to gain and phase margins of -6 db - 18 db and ±52 degrees 
respectively. This is a considerable improvement ever the 
original design. The factor that changed the design was the 
upper diagonal feedback term which provides a cancelling 
factor for the cross-coupling term b ^ = 50. This can be 
seen by looking at the system matrix equation 6.11. 



A = A - BF 

C «'» % 



where equation 6.12 gives A 



A 




~ ^ 12. 12 ^2Z 



-1-f 



Z2 



( 6 . 11 ) 



( 6 . 12 ) 
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Figure 6.8 Singular Value Plot Case Three. 
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If the upper diagonal term of the matrix can be driven to 
zero the system will be decoupled into a diagonal system. 
To place the poles cf this system at - 2,-2 the diagonal 
feedback gains must he one. The optimizer feedback gains 
for this case were f /; =1.00 006 .99998 and, most impor- 

tantly/ f ,£ =-5 1 . 834 0 1 . Using these values in the A matrix 
A = -2.00006 1.83501 

0.0 -1.99998 

The value in the upper right position in the system matrix 
has been lowered considerably from the value of near 50 that 
appears in this position in the low singular value cases, 
lowering this system gain value decreases the cross-coupling 
perturbation effects cn the system. This can be more graphi- 
cally demonstrated by figures 6.9 and 6.10. In figure 6.9 
the transfer function shows a high gain of approximately 35 
db and a bardwidth of 50 rad/sec. In figure 6.10 this gain 
has been reduced to 6.0 db with a bandwidth of almost 2 
rad/sec. Two things are indicated by the figures, one, the 
open-loop Ecde plot of the transfer function of the cross- 
coupled channel can be used to indicate the robustness 
problem as evidenced by the high gain and bandwidth relative 
to the other feedback gains and, two, the mechanism used by 
the pole placement and robustness design procedure to 
recover robustness is to reduce the relative gain and asso- 
ciated bandwidth within the affected channel. Figures 6.11 
and 6.12 which are for the input one to output one channel 
are approximately the same as are the Bode plots for the 
input two to output two channel which are not shown. These 
figures indicate that no problem exists in the diagonal or 
direct coupling terms. The pole-zero diagram for optimiza- 
tion run case three is included in figure 6.13. The only 
significant change in this plot as compared to figure 6.7 is 
the movement of the zero in the off diagonal pole -zero plot. 
The zero is seen to shift to one of the pole locations 
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Figure 6.9 Open-loop Transfer Function 2-1 for Baseline. 
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Figure 6.10 Transfer Function 2-1 Optimized! 3 Var). 
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figure 6. 11 transfer Function 1-1 Baseline. 
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Figure 6.12 Transfer Function 1-1 Optimized (3 variables). 
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Figure 6.13 Closed- loop Pole-Zero Plot Case 3. 
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during the optimization. This zero shift has the effect of 
smoothing the freguency response curve in the vicinity of 
the freguency of the minimum singular value providing a more 
uniform gain distribution around this point. Table 1 shows 
comparison results for the feedback gains for several cases 
of this fcasic problem. 

The logical extension of the problem to the general case 
is to allow all four feedback gains to become design vari- 



T ABLE 1 

Comparative Results Simple Problem 





Feedback Gains 




f ll 


f 12 


f 21 


f 22 


f ^2 variable f 2 l = ® 


1.00006 


-51.B3401 


o 

« 

o 


0.9998 


f 21 variable f 12 =0 


0.25254 


0.0 


-0.01071 


1.52332 


f^^ variable f^ 2 = ~50 


.99960 


-50.0 


0.00025 


0.99038 


All f's variable 


1.11936 


-60.4718 


0.00280 


0.75478 



ables. The pole placement and robustness routine can then 
use full freedom in choosing all four of these feedback 
gains to compensate for any cross-coupling effects within 
the system. Based on the previous analysis the two diagonal 
gains would be expected to approach 1 while the upper off 
diagonal gain moves tc -50 and the lower one moves to 0 . The 
flexibility of the ACS program was reguired for the four 
design variable study. Several runs were made with various 

combinations of starting parameters for the feedback gains 
and optinizer routines before a good design for the case 
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incorporating all four feedback gains as variable was 
obtained. While this points out one of the limitations of 
optimization design routines the program was able to develop 
an improved design ever the baseline while employing four 
design variables. In the design produced for this case the 
singular value level was placed above 0.6 for a gain margin 
of -4 db to 9 db and phase margin of ± 35 degrees. In this 
formulation the optimizer was able to place the poles and 
meet the design singular value level. The feedback gains 
produced by the optimizer were: 

1.11936 -60.4718 

0.0028 0.75478 

which are approaching the analytic design gains. After 
obtaining the feedback gains the OPTSYS program was used to 
obtain the necessary data to 'do a closed-loop pole-zerc map. 
This plot is shown in figure 6.14 These plots indicate a 
similar pole-zero location to that found in the previous 
three design variable problem. Again the gain in the 
affected channel has been reduced to compensate for the high 
cross-coupling perturbation within the system. 

The design studies presented to this point have been 
based on breaking the system loop at the input as shewn in 
figure 6.15. In multivariable theory the location of the 
break in the loop changes the return difference for the 
system and the transfer function formulation. In the figure 
number 1 depicts a system with an input loop break point 
while number 2 is an output loop break point for output 
return difference determination. The return difference func- 
tion for the point 1 is written as I+Fg while the return 
difference fer point 2 is I+GF. The baseline system not only 
has low singular values of the input return difference 
matrix, the lowest being 0.0706, but also has low singular 
values for the output return difference matrix of point 2 of 
about the same order. To demonstrate the versatility of the 
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Figure 6.14 Pole-Zero Plot Simple Problem (4 fariables). 
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Figure 6.15 System Block Diagram. 

pole placement and robustness routine a design was made in 
which the output singular values were specified instead of 
the input values. The pole placement and robustness routine 
produced a design with gains that provided a highly robust 
system on the output side with gain and phase margin cf -6 
db to and phase margin of ±60 degrees. The design was not 
robust on the input side. Thus, designing for robustness at 
one point in the system does not necessarily give robustness 
at all points within the system. 

Cne final case that should be discussed is that of 
setting both input and output robustness criteria at the 
same time. Excellent results were obtained for this case. 
The design routine placed the poles at -1.97±.009j with 
feedback gains of 

0.85793 -45. 91757 

0.00425 0.87137 

The input singular values were raised to a level of 0.74984 
or -4.5 db to 12 db gain margin and ±43 degrees phase 
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margin. The output values were raised to above 0.822 which 
corresponds to -5 to 15 db gain margin and ±49 degrees phase 
margin. 

Tc summarize, it can be stated that the robustness 
problem for this system exists in the upper cro ss- coupling 
channel (input two,output one). The lack of robustness can 
be discovered in two ways. The first method is to plot the 
open-loop Bede plots cf each element of the transfer matrix 
and lock for extremely high gains and bandwidths relative to 
the ether transfer functions. The second method examines 
the singular values of the return difference matrix for 
magnitude. Low singular values correspond to low robust- 
ness. The pole placement and robustness design routine can 
increase robustness by modifying feedback gains tc reduce 
the effect of cross-coupli ng within the system. Observing 
the gain modification made by the pole placement and robust- 
ness rcutine the critical channel within the system that 
affects the robustness may be determined from the Bede 
plots. The pole placement and robustness routine feedback 
gain changes also cause zero shifts during the robustness 
recovery. The gain on the open loop Bode plot for the 
affected crcss-coupling channel is adjusted and the clcsed- 
loop zeros, as seen cn the pole-zero diagram, are shifted. 
This zerc shift is in a direction which will combine with 
system pcles to smooth the freguency response diagram in the 
vicinity of the minimum singular value. 
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VII. A HEIICOPTEB STA BILI TY PROBLEM 



This chapter will' deal with a more practical application 
of the numerical optimization proyram. In this problem the 
combined pole placement, robustness design procedure will be 
applied to the linear lateral dynamic channels of a CH-47 
helicopter. The model is a highly coupled two-input two- 
cutput system that has been studied for its basic robustness 
characteristics £Hef. 16]. The usual procedure for design 
of highly coupled systems is to obtain a diagonally dominant 
closed-lcop system. This diagonal system will be stable but 
not rcbust to cross-feed parameters. Sandell,et al, produced 
three designs. Two cf these designs, while meeting basic 
performance criteria, had poor robustness. The third design 
was a relatively good design. The numerical optimization 
technique developed in the thesis was applied to the two 
poor designs and shown to provide substantial improvement in 
robustness. 

The systems were designed to satisfy specifications to 
step input response and stability margins as stated in mili- 
tary specifications. Specific design parameters for each of 
the three designs presented were not available. The designs 
were all considered to meet the performance specification 
criteria and stability margin requirements. It was shown 
that two of these designs were extremely sensitive to model 
errors. The classic Nyquist techniques did not predict this 
sensitivity. The singular value analysis did indicate sensi- 
tivity problems. 

The model is that of a CH-47B helicopter lateral dynamic 
system in hover. The dynamic model of the system is 
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X = A X + 3 U 
— /> * — 



(7.1) 



X =(v,p,r,y) 



(7.2) 



u 



“ (Sir Sc) 



(7.3) 



wher e 



A = 



2.27 


-1.42 


-0. 1 5 


31.99 


0.01 


-0.7 


-0. 07 


0 


0.04 


-0.05 


-0. 5 


0 


0 


1 


0. 11 


0 



and where 



B = 



0. 12 
0.04 
.34 
0 



0.9 5 
-8.37 
0 . 02 
0 



with full state available for feedback. Table 3 is a summary 
of parameters. The system is not open-loop stable. 

Three control laws are formulated to satisfy the desired 
performance specifications. Equation 7.4 is the basic 
control law. 



u = -F -x + h-. 

• L t T 

i= 1, 2 or 3 

where the following values of F and h were are; 
-1.72 -23.5 70.6 595. 

0. 024 -2 .71 0.368 -7.99 



(7.4) 



= 



= 



- 



0. 198 154.0 18.3 142.0 

-0.01 -1.592 -0. 139 -1.47 



0 

0 



0 

-4 



25.5 

0 



0 

-27 
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TABLE 2 

CH-46 Helicopter Parameter Definitions 



Variable 


Units 


Description 


V 


ft/sec 


Vehicle body y-axis earth relative velocity component 


P 


rad/sec 


Roll rate 


r 


rad/sec 


Yaw rate 


* 


rad 


Roll attitude angle 


6 b 


in 


Yaw rate rotor deflection control 


6„ 


in 


Roll rate rotor deflection control 


c 








The is a step input command that must track. Figure 

7.1 is a diagram of the control structure for the feedback 
control laws. Figure 7.2 is a detailed layout of the 
system. All three cf the designs have negative real eigen- 
values and provide stable overdamped responses. All three of 
the above designs meet the desired design specifications. 

The loop Bode plots, £Bef. 16], indicate all three 
designs to have good stability margins when considered a 
loop at a time. Since these designs were all full state LQ 
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Figure 7.1 Feedback Control Structure. 
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Figure 7.2 System Diagram. 
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designs with a diagonal control weighting matrix they should 
posses at least -6 db to infinite gain margin and 60 degrees 
of phase margin. fchen the singular values of the three 
designs are computed a robustness problem is indicated by 
low singular values cf two cf the designs. Figure 7.3 pres- 
ents singular value plots of all three designs. This plot 
shows that designs 1 and 2 both have very low minimum 
singular values for the return difference matrix. Design 1 
goes as low as -20 db near 10 rad/sec in frequency while 
design 2 is down to -34 db at frequencies up to 1000 
rad/sec. Design 3 is a good design with singular values that 
remain above one throughout the frequency range of interest. 
Using the universal gain and phase diagram as discussed in 
Chapter 6 this equates to a gain margin of about -6 db to 
infinity and a phase margin of 60 degrees as expected. 

The system stability in design 1 may be affected by 
perturbations occurring in the actuators as shown in figure 
7.4. If the output axis coupling from £c spills into in 
the frequency range from 0.5 to 50.0 rad/sec with a magni- 
tude cf 0.12 and a phase of 60 degrees then the system can 
become unstable. This could be caused by nonlinear terms, 
worn parts, or system saturation. 

The second design as shown in figure 7.5 may have a 
stability derivative variation between Sq and p. If this 
derivative varies from about 0.04 to -0.96 the system can 
become unstable. The center of gravity location, trim of the 
aircraft and rotor coupling can all affect this derivative. 

The pole placement and robustness design technique 
developed for this thesis research was applied to designs 1 
and 2 to obtain an improved robustness for these designs. 
Design 1 will be considered first. In this design a 

cross-feed perturbation through the actuator can produce 
instability as indicated by the low singular values at lower 
frequencies. To study this problem with the pole placement 
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figure 7.3 Alphatech Design Singular Value Plots. 
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Figure 7.4 Design One Perturbation In pot. 
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Figure 7.5 resign Two Perturbation Input. 
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and robustness design method it is assumed that the eigenva- 
lues (pole locations) of the system as developed in 
[Ref. 16]. are the required poles for the performance 
criteria. Once the pole locations are set, robustness 
criteria must be selected. Prom the universal gain and phase 
margin curve discussed earlier two choices of singular value 
levels were made for this problem. The first singular value 
level chosen was 0.6. This corresponds to a gain margin of 
-4.0 db to 8 db and a phase margin of about 35 degrees. The 
second value chosen was 1.0. With corresponding gain margin 
of -6 db toco and a phase margin of 60 degrees which are 
the characteristics of a LQ regulator design with diagonal 
weighting matrices. 

For singular value level 0.6 the pole placement and 
robustness design routine places the poles as shown in table 
2 The slight differences in these pole locations appear to 
have insignificant effect on the performance as shown in the 
response curves of the system. The feedback gain adjustment 
moves the minimum singular value from about 0.1 1 with very 
poor phase' and gain margins to a level of 0.66. This is 
above the desired values of gain and phase. The improvement 
in robustness came from modifying the feedback gains in 
channel Ss . 3y greatly reducing the gains in channel Ss the 
optimizer minimizes the influence of the cross-coupling from 
channel Se. . In this way a much larger spill over of channel 
Sc may be tolerated through the actuator before the system 
will become unstable. 

The primary mechanism of robustness improvement in this 
problem was a reduction in the gain levels of the affected 
channel. The feedback gains presented in table 4 show the 
modification of these gains from those utilized in 
[Ref. 16 ]. While all the gains are modified the f gain 
undergoes a significantly larger change than the ether 
design 1 case 1 feedback gains. Looking at the open-loop 
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TABLE 3 

Design One Pole Placement 



Pole 


Desired Location 


Actual Location 


1 


- 24.7977 


-24 . 7893 


2 


-11.3635 


-12.0083 


3 


-10.3288 


-10.7752 


4 


- 2.1005 


- 2.1181 



transfer functions cf this optimized problem and comparing 
them, with a non- optimized set of transfer functions for S& 
to as shown in figure 7.6 and 7.7 show very little 

change in the system gain, both are around 42.0 db, and a 
bandwidth of about 20 rad/sec. There is a significant 
change in the phase diagram which is caused by the zero 
location shift. 

The transfer function for o>c to Sc. depicted a gain 
increase of about 3 db for the optimized design while the 
phase remained similar for both transfer functions. In 
transfer function Sc. to S& some important aspects of the 
problem are observed. The Bode diagram of the open-loop 
transfer function of Sc to S& clearly indicates the 
cross-coupling problems and the pole placement and robust- 
ness design routine's mechanism of optimizing the system 
gains to increase robustness. The gain is reduced from 93 
db to about 78 db and the bandwidth is reduced from above 
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Figure 7.6 Transfer Function £ q~£q Design One, Nonoptiaized. 
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Figure 7.7 
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100 rad/sec to about 30 rad/sec. The reduction of bandwidth 
and gain in the loop Sc. tc Sb yields an increased toler- 
ance to perturbation. The transfer function for Sc- to 

<Ss is shown in figures 7.8 and 7.9. In the transfer func- 
tions fox S& to £c- the bandwidth is slightly increased 
from 0.6 to 0.8 rad/sec and the gain actually increased from 
9 db to 14 db. The big change in the overall system, 
however, is in the transfer function from S<~ to Sa . This 
is the channel that the destabilizing perturbation enters 
and by greatly reducing the gain and bandwidth in this 
channel through a change in feedback gain, the optimizer 
routine has brought the entire system gains to more balanced 
conditions and recovered a highly robust design. 

The gain changes associated with the robustness improve- 
ment cause the zeros ox the various closed-loop pole-zero 
diagram tc move. A comparison of the eight pcle-zero 
diagrams is shown in figures 7.10 to 7.13. The significant 
feature of these pcle-zero diagrams is the shift of the 
zeros of the optimized design in a direction that attempts 
to equalize or balance the frequency response for frequen- 
cies in the vicinity of the minimum singular values. The 
pole-zero diagram of Sb to v will be discussed as an example 
of this effect. In 7.10 the nonoptimized zeros are located 
about -2 ± 4. 5 j and -19.8. When the pole placement and 

robustness routine has completed the feedback gain modifica- 
tion these zeros have shifted near -11 and -6±3 j . The 
effect of these zero shifts is to combine with the pcle 
locations tc equalize the frequency response as depicted in 
figures 7.14 and 7.15. Zero shifts for the remainder of the 
transfer functions provide similar results in the ether 
channels. By moving toward the frequencies associated with 
the minimum singular values the zeros have balanced the 
overall frequency response of the system in each channel. 
While the channel gain modification is the primary mechanism 
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Figure 7.8 
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Figure 7.9 Transfer Function S c ~ S3 Design One, Optiaized. 
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figure 7.10 Pole-Zero Plots for Design 1. 



98 



NON OPT 



OPT 




8 




V 





©— < — 3—1 — ■ — I — • — I 



-Zo -to 



JO £0 



8 



c 



* p 



5 C " p 



Figure 7.11 Pole-Zero Plots for Design 1 (cont.) . 
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Figure 7.12 Pole-Zero Plots for Design 1 (cont.) . 
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Figure 7.13 Pole-Zero Plots for Design 1 (cont.). 
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Figure 7.14 The S 3 to v Freguency Response, Nonopt.. 
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Frequency Response, ILptiaized. 
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FREQUENCY - RAD/SEC 



for robustness recovery the zero shift associated with the 
feedback gain changes is directly related to the overall 
freguency response of the system. 

In design 1 case 2 the minimum singular value level was 
increased from 0.6 tc 1.0. The pole placement and robust- 
ness design procedure also gave a better design for case 2 



T1BLE 4 

Helicopter Problem Feedback Gains 



EE SIGN 


GAIN VALUES £ll *12 *13 Jl4 
.21 22 I 23 I 24. 


One 

Case 1 


-14.77726 2.15858 77.96629 -32.91595 

-0.00567 -2.55646 0.39039 -15.04805 


One • 
Case 2 


-6.36249 -1.53746 72.85013 74.20387 

0.00813 -2.60073 0.65395 -13.47128 


Two 


-0,96660 -4.54597 0.75436 11.96808 

0.00916 -3.03069 -0.36391 -0.34991 - 



by modifying the feedback gains as shown in table 4. The 
change in the transfer functions for Sc to Sg {figure 
7.16 and 7.17) is similar to the change seen in figures7.8 
and 7.9. The gain was reduced from above 90 db to about 75 
db and the bandwidth cut from above 100 rad/sec to about 25 
rad/sec. The added requirement of increased robustness did 
not result in a significant change in performance. 
[Ref. 16 ]. 
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Figure 7.16 



Transfer Function 




Design 1 Case2, Ncnopt.. 
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Figure 7.17 Transfer Function $ c - Design 1 Case 2, Opt. 
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Design two which is the design that is subject to param- 
eter uncertainty in the to p channel was also studied 

using the pole placement and robustness design routine. In 
this design the singular values were very low until frequen- 
cies above 100 rad/sec were reached. To improve the design 
it was assumed that the pole locations that corresponded to 
[Ref. 16] were the required pole locations and that a 
robustness singular value level of 0.6 would provide 
adequate gain and phase margin for the design. The pole 

placement and robustness routine adjusted the gains in this 

problem until a minimum singular value of 0.6 was obtained. 
During this adjustment the gains in channel were consid- 
erably reduced to offset the cross-coupling between the two 
channels. Again, in this design as before the channel 
cross-coupling between and s* shows a marked change 

in bandwidth and gair from above 100 to 35 rad/sec and about 
78 db to 22 db respectively. Figures 7.18 and 7. 19 show the 
transfer function plcts for this term. Figure 7.20 shows 
the singular value improvement. In figure 7.21 the time 

response is plotted. As can be seen in the plot the 

improvement in robustness for this problem results in very 
sluggish response and degraded performance. 

To summarize for this problem, a given performance level 
has been chosen in terms of pole locations. The level of 
robustness has been set for a desired gain and phase margin 
based on the universal gain and phase margin curve. The pcle 
placement and robustness routine has been able to improve 
the robustness level by changing the feedback gains that 
affect the channel cross-coupling. This robustness 
recovery is affected by modification of the system feedback 
gains in such a manner that cross coupling gains are reduced 
so that small cr oss-ccupling perturbations do not drive the 
system into instability. The open-loop transfer function 
plots have teen used to indicate how this mechanism operates 
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Figure 7.18 Transfer Function S c ~ Sa Design Two, Nonopt.. 
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Figure 7.19 



Transfer Function S c - St Design Two, 



Optimized. 
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SINGULAR VALUE 
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Figure 7.20 Singular Value Plot Design Two. 
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Figure 7.21 Tine Response Design Two. 



Ill 



and have teen shewn tc be an alternative indicator of chan- 
nels that may be affected by cross-feed perturbations. Ihe 
pole-zero diagrams of the closed-loop transfer functions of 
the transfer matrix further indicate that zero movement is 
in a direction that equalizes the gain level of the 
frequency response curves in the vicinity of the lowest 
singular values providing a more balanced system response. 
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VIII. SiaPIE O BSER VER 



The pcle placement and robustness design procedure can 
also be used for robustness recovery in observer design. 




Figure 8.1 Simple Observer. 



Given an observer as represented in figure 8.1 it has been 
shown that the system differential eguation may be written 



[uation 


8. 


1 












• 

X = 


p 




-BF* 




+ 


V 


r 


(8.1) 


A 

X 


[KC 


Ac 




k 




B 






*c 


= A 


- 


KC 












■ X 


is 


the 


sta te 


and 


X 


is the observer variable. 


and the 



where 

solution may be separated to independent solutions for the 
feedback gains and the observer gains. This separation 
allows the use of the feedback gains to set the pole 



113 



locations and then a second optimization run with set feed- 
hack gains so that the observer gains may be computed to 
adjust the system robustness level. The observer pole loca- 
tion could also he placed using the observer gains, K. In 
the currently implementation of the pole placement and 
robustness design routine the observer poles are simply 
restricted to areas cf the left half plane. 

In this chapter a simple stable observer system will be 
analyzed based on [Bef. 17]. Given the system, equation 
8 . 2 , 



x = 



Fo r 


X + 


m m 

0 


u + 


35 


-3 -4 




1 




-6 1 
J 



( 8 . 2 ) 



J s l 2 1 J X + rj 

where ECf)=E( n )=0 and E ) 7 (r) ) =E ( t ) 1 (r) ) = t-f) 
the feedback law of equation 8.3 



w ith 



u=-(5C 10) x + SOr 



(8.3) 



An analysis has been done to compare results of the numer- 
ical optimization procedure with results presented in 
[Ref. 17 ]. The results for an optimal regulator design 
using quadratic cost criteria as stated in equation 8.4 are 
given in table 5 

J = J (x T H T Hx+ u 2 } dt (8.4) 

with H = 4 y? ( ^35 1) 

Figure 8.2 shows a Nycuist plot of the full state regulator, 
the optimal filter and a fast filter. The full state design 
had poles at s= -7. 0±j 2.0 and feedback gains of 50 and 10. 
The optimal filter as shown in table 5 had poles of -7. 0 ± j 
2.0 with gain and phase margins of -6.75 db and ±15 degrees. 
The optimum filter gains were 30 and -50. Using a faster 
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Filter with 

Poles at 

-22+17 .86 j 



Optimal Filter 




Full State 



Figure 8.2 Nyguist Plot. 

filter also gives peer gain and phase margins. The gain 
margins are on the order of -.98 while the phase margin is 
less than 10 degrees. The bandwidth also increased from 12 
to 40 rad/sec. A recovery procedure based on a modification 
to the process noise matrix [Ref. 17] may be applied to the 
problem. This procedure can recover a large amount of the 
robustness that was lest with the observer addition. Figure 
8.3 shows data obtained for several trials of tha fictitious 
noise procedure. The gain, phase margin and ether 
parameters may be found in table 5. 



TABLE 5 

Observer Parameter Data 





FILTER 


GAIN 


PHASE 


ERROR 


STATE 




rum 




ROLLS 


KARGIN 


KARGIN 


COVARIANCE 


covariance 




UIX 






db 


deq 


l ((*■ 




Ef*« 7 ) 






Optical 


.J i2) 






$7 


-163 


221 


-613 


30 


LOG Design 


- 6.75 


15 


-163 


277 


•613 


2070 


-50 




Fast filter Adjust- 








6280 


-12200 


no 


-613 


720 


ment Procedure 


- 22 * 77 . 9 ) 


- .98 


10 


-n:oo 


23800 


-613 


8520 


-1400 


Fictitious »olse 


•1.3 






107 


-184 


236 


-613 


2i.i 


Adjustment Procedure 




- 7.73 


19 












q* • 100 


-11.1 






•184 


319 


-613 


1810 


00.1 




-2.9 






163 


-301 


263 


-«u 


20.4 


q‘ • 500 


•24 


-10.9 


33 


-301 


564 


•613 


uoo 


-17.7 


q* • IQ* 


•2.5 






204 


•MS 


285 


-613 


16.7 


1 1 w 


-33 


-13.9 


42 


-385 


743 


-613 


1360 


•1.9 




-2.1 


-37 


74, 


290 


-570 


317 


-613 


6.9 




-100 






-570 


mo 


-613 


1200 


34.6 



POPLAR CASE I 1.00018 

-2.02742 



GM=-5.5 dB 
PM=117 deg 



POPLAR CASE II 19.22148 

32.59476 



GM=-19 dB 
PM=58 deg 



The singular values of the state feedback and optimal 
observer systems were computed for comparison / ith results 
produced using the pole placement and robustness design 
recovery procedure. Figure 8.4 indicates a loss in robust- 
ness represented by the observer singular values. lower 
singular values are less robust. For this 5 ingle- input 
single-output observer the Nyguist diagram will be used to 
define gain and phase margins. Since a singular value of 1.0 
is indicative of a linear quadratic level of gain and phase 
margin, i.e. GM=-6 db, and PM=±60 degrees, this was chosen 
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figure 8.3 Hyguist Plot for Robustness Recovery. 

as the design level. The pole placement and robustness 
routine vas used to recover robustness while setting the 
pole placement at the state feedback pole locations of 
-7±j2. The minimum input singular value level was set at 
1.0. The pole placement and robustness routine was also set 
to place the observer poles anywhere between -100 and -2 
that would provide robust design. The plant poles were 
placed at -7.05 ±j 1.82. The feedback gains for this run 
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SINGULAR VALUES 




Figure 8.4 Siagular Values of Observer Systea. 
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were 50 and 10.09. The optimizer output produced a filter 
with a gains of 1.00018 and -2.02742. The observer pcle 
locations were -1.99 ± j .001, very near the plant zero 
location. A significant change was produced in singular 
values. Figure 8.5 shows all singular value curves plotted 
together. The optimizer solution for this problem is well 
above the optimal filter curve at low frequency. Figure 8.6 
shows the Nyquist plot of the optimizer developed design. 
The system has a gain margin of -6 db and a phase margin of 
117 degrees. The most significant differences between the 
two designs being that the observer poles are close to the 
plant zero locations and the filter gains are much lower for 
the optimizer solution. 

Using the OPTSYS program with the pole placement and 
robustness routine computed gains as the design parameters 
the data for the observer filter was computed. The error 
covariance matrix was found to be 
306.5 -456.5 

-456.5 7C0.1 

These values compare favorably with the trends established 
in table 5. The last comparison of the pole placement and 
robustness design recovery procedure was an analysis of the 
time response curves. Figure 8.7 shows the comparison plot. 
The design obtained using the pole placement and robustness 
routine did not degrade system performance. 

Cne additional analysis was conducted that set the 
desired design parameters slightly differently. In this run 
the pcJe placement and robustness routine was set to place 
the observer poles between -10 and -100. The pole placement 
and robustness routine was unable to totally satisfy this 
requirement. It violated one of the constraints and moved 
the larger observer pole to -2.25 which is near the optimum 
pole location and also the high q 2 values of the fictitious 
noise procedure. The smaller pole was moved out to -70.8 
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Figure 8.5 Singular falue Comparison Plot. 
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Figure 8.6 Nyguist for Computed Robustness Recovery. 



which corresponds to the same level of pole movement found 
when g 2 becomes large. The singular value levels were raised 
significantly as shown in figure 8.5. The state and error 
covariance matrices were of the same order as the matrices 
found in table 5. This design has a phase margin cf 58 
degrees. The gain margin was about -19 db. Even though seme 
of the constraint conditions on the design were not met the 
design demonstrates excellent robustness recover/ . 
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Figure 8.7 Time Response Plot for Simple Observer. 
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The pcle placement and robustness design routine has 
produced a robust design for this observer based system. 
This design is obtained by using a numerical optimization to 
directly manipulate the feedback and filter gains. 
Modification of the IQ functional equation as done in the 
fictitious noise procedure is not required. The pole place- 
ment and robustness routine solution for this problem has 
lower gains than those found using the fictitious noise 
adjustment. The routine provides a good, direct methodology 
for selecting the feedback and filter gains for a robust 
observer design with excellent performance. 
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IX. E03DS T 03 SE 57 SR DESIGN 



This chapter will be devoted to a short discussion of 
the robustness recovery of a fourth order observer based 
controller. The problem is a helicopter problem, [Ref. 6]. 



In this 


case 


the helicopter model 


is that 


of 


the 


longitu- 


dinal control 


loop 


of a CH- 


47. The 


nominal 


model 


is 


taken to 


be the system 


of equations 


9.1 and 


9. 


2 for 


an aircraft speed 


of forty knots 




















• 
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’-0. C2 


0.005 2.4 


-32. 
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0. 14 
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u 


(9.1) 




-0. 14 


0.44 


-1.3 


-30. 






. 36 
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0 . 


0.018 -1.6 
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.35 


.009 








0 . 


0 . 
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■0 
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0 1 0 
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(9.2) 




0 0 0 


57.3 


















In this 


problem the contro 


Her 


is 


formulated as 


in figure 



9.1 which leads to an open-loop 
form of equation 9.3 



transfer function of the 



K (S) G (s) =F ( SI-A+ EI+KC) ” 1 KC (sI-A) -IB 



(9.3) 



The pole placement and robustness recovery procedure was 
applied to this problem. 

First, the standard full state feedback design was 
carried out using the Naval Postgraduate version of OPTSYS. 
This design produced excellent singular value output for the 
return difference as shown in figure 9.2. The lowest 
singular value being essentially 1, corresponding to a LQ 
design with -6 db to cd gain margin and 60 degrees of phase 
margin. The time response of the system was good as shown in 
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Figure 9. 1 Observer Based Controller. 

figure 9.3, reaching steady-state in about four seconds with 
only a slight overshoot. Since the characteristics for this 
system are acceptable no further design iteration was 
carried out. The full state feedback became the baseline 
design. Assuming that full state feedback was not available 
and only, two measurements could be produced an observer was 
developed to control the two measured outputs, vertical 
velocity and pitch attitude. Using the measurement matrix of 
equation 9.2, OPTSYS was used to develop an optimal observer 
system for this problem. The singular values of the return 
difference matrix and the time response were plotted for 
comparison to the full state design. These are shown in 
figures 9.2 and 9.3 respectively . Note that the singular 
value of the return difference matrix is as low as 0.16 at 4 
rad/sec. This equates to a gain margin of -1 db to 1 dt and 
less than 10 degrees phase margin. The time response is 
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Figure 9.2 Singular Value Plot of Observer Results. 



126 




9X0 



900 



92*0- 



09 0- 

2X 



94*0- 



00*t- 9X1- 



r— 

9*0 



— r 
CO 



— , — 

ro 

CX 



— r- 
0*0 



— I— 

x*o- 



i 

9*0- 



*0 



2*0 



t*0- 



I '' ' I ' ' f . m. . j r 1 i I 1 

09*0 92*0 02 0 91*0 01*0 90*0 00*0 90 0- 

YX 



Figure 9.3 Time Response for Observer System. 
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TIME - SEC 



plotted in figure 9.3. The steady-state was reached for the 
velocity after a slight overshoot in about four seconds. 

As the final step in the analysis the pole placement and 
robustness routine was employed to recover robustness of the 
observer based system. The pole placement and robustness 
routine was first used to set the poles at approximately the 
same location as the poles of the LQ regulator. These were 
assumed to be the desired pole locations. The pole placement 
and robustness routine was then used to vary the filter 
gains, K, until the desired level of robustness was reached. 
The desired singular value level chosen was 1.0 which corre- 
sponds to -6 db toco db gain margin and 60 degrees of phase 
margin. Figure 9.2 shows that the pole placement and robust- 
ness procedure failed to totally recover the robustness 
level to 1.0. The minimum singular value reaca ed was only 
0.936. This singular value equates to a gain margin of -5.0 
db to a» and a phase margin of 55 degrees. While this is 
slightly less than the design objective, it is Ear superior 
to the optimal observer design discussed in the previous 
paragraphs. This recovery was made by making an optimiza- 
tion run, finding the filter gains with near zero values and 
freezing these values at zero. This reduced the number of 
design variables the pole placement and robustness routine 
was required to manipulate in a second optimization run and 
gave a higher robustness solution. 

This analysis and the second order observer analysis 
presented earlier clearly indicate that it is possible to 
use the pcle placement and robustness procedure and separa- 
tion principle to develop a robustness recovery procedure. 
The pole placement and robustness routine provides robust- 
ness recovery by direct modification of the feedback and 
filter gains. This procedure requires no modification to LQ 
cost functionals or ether parameters as done in the ficti- 
tious noise adjustment method commonly used for robustness 



128 



recovery. By providing direct gain adjustment the pole 
placement and robustness procedure results in a practical 
design with relatively low observer gains and good perform- 
ance. The procedure is simple and straight forward with the 
only difficulty being a requirement to sometimes modify 
initial starting values or optimizer codes to force the 
solution toward the desired point. 
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X. CONCLUSIONS 



An effective method of robustness multivariable control 
design utilizing a numerical optimization based algorithm 
has been developed. The pole placement and robustness design 
routine coupled with the Automated Design Synthesis program 
provides the designer an excellent tool with which tc attack 
the robust design problem. 

The pole placement and robustness design routine has 
demonstrated the capability of providing designs that solve 
the problems caused by cr oss-coupling perturbations which 
reduce robustness in multivariable systems. This design 
improvement is accomplished by modifying the system feedback 
gains in such a manner that the gain in channels that are 
affected by cross-coupling perturbations is equalized with 
other system gains to reduce this cross-coupling effect. The 
gain changes are accompanied by zero shifts which also 
influence the gain distribution and frequency response of 
the system. 

Perturbation problems in multivariable systems have been 
shown to be detectable by singular value analysis and by 
using the Bode magnitude diagram of the open- loop transfer 
functions of the system. In the open loop transfer function 
large differentials in Bode gains and bandwidths are indica- 
tive of problem areas for cross-coupling perturbations. 
Robustness is obtained by the pole placement and robustness 
design program by modifying those gains and bandwidths asso- 
ciated with the cross-coup ling perturbations thus reducing 
the amount of energy coupled from the perturbation into 
other channels. An associated zero shift has been observed 
when these gain modifications take place. This zero shift is 
in the direction of poles that are located in the vicinity 
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cf the frequency of the minimum singular value and tends to 
equalize the frequency response curve gains in this region. 

The use of numerical optimization to recover robustness 
in observer based designs was demonstrated. The pole place- 
ment and robustness routine was applied to problems 
previously solved using the fictitious noise procedure for 
robustness recovery. The direct manipulation of feedback and 
filter gains by the pole placement and robustness routine 
provided a highly robust design with relatively lew filter 
gains. The problem cf robustness recovery in filter- 
observer designs has been solved in a straight forward and 
highly practical 



manner. 
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